First order wetting of rough substrates and quantum unbinding 
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Abstract 

Replica and functional renormalization group methods show that, with short 
range substrate forces or in strong fluctuation regimes, wetting of a self-affine 
rough wall in 2D turns first-order as soon as the wall roughness exponent ex- 
ceeds the anisotropy index of bulk interface fluctuations. Different thresholds 
apply with long range forces in mean field regimes. For bond-disordered bulk, 
fixed point stability suggests similar results, which ultimately rely on basic 
properties of quantum bound states with asymptotically power-law repulsive 
potentials. 
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Wetting transitions occur when, e.g., an interface separating two coexisting phases un- 
binds from an attractive substrate. In the recent literature, much space has been devoted 
to the effects of different types of disorder on the nature and universality of such transitions 
This is motivated both by the presence of impurities in actual experiments, and by 
the expectation that disorder modifies critical behavior. Many studies concentrated on im- 
purities in the bulk, or on the surface of a smooth substrate. Another type of disorder is that 
due to the roughness of the wall delimiting the substrate. This rather frequent geometrical 
disorder was discussed especially in connection with measurements of nitrogen adsorption 
on flash deposited silver 

Substrate roughness describable by self-affine geometry is often realized and most in- 
teresting, from both a fundamental, and a physical point of view. Indeed, a wall whose 
average transverse fluctuation, Wl, increases as a power of the longitudinal sample size, L 
(Wl oc L^ w ,0 < ( w < 1), has a random geometry characterized globally by a single roughness 
exponent, ( w . Moreover, no matter whether their fluctuations are controlled by temperature 
or by disorder, bulk interfaces behave as self-affine objects, with appropriate exponents (q 
]Tjrj]. Thus, one may expect a direct competition between wall and interface roughnesses to 
take place in wetting phenomena, when substrates of this kind are considered. 

In the present Letter we show exactly in 2D that, as soon as the wall wins (i.e. for ( w > (q 
with short range substrate potentials), the above competition is resolved in an unusual, 
drastic change of the wetting transition, from continuous to first-order. We determine exact 
roughness thresholds for first-order wetting also in other regimes with the substrate exerting 
long range forces (e.g., van der Waals) on the interface and discuss, at perturbative level, 
cases with bond-disorder in the bulk. A change from continuous to discontinuous wetting 
is in fact a quite surprising and unexpected phenomenon, especially in 2D. As a rule, with 
short-range forces, first-order wetting never occurs in 2D, except in special and ad hoc 
limit situations [[J]. Disorder generally acts in the opposite sense of turning into second- 
order discontinuous transitions |TJ . In all cases the mechanism leading to first-order wetting 
at large ( w can be traced back to some basic properties of quantum bound states of a particle 
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in a potential |IQ |. Such properties turn out to have very far reaching consequences for our 



understanding of wetting phenomena. 

In 2D, a self-affine wall can be described by a random function h w (x) (Fig. 1), with 
probability distribution j7| 

P w [h w \ oc exp [- J dx^^K/dx?) 2 ] (1) 



such that \h w (x) — h w (x')\ oc \x — x'\^ w , with ( w = (3 — 1/2 and the overbar indicating 



quenched average with respect to (1) [11|. In the presence of bulk disorder, the interface 
Hamiltonian can be put in the form 

(3H[h, h w , V] = J dx[^(dh/dx) 2 + U(h(x) - h w (x)) + 

+V(x,h(x))} (2) 



where V is a Gaussian random potential (V = 0, V(x, h)V(x', h') = A5(x — x')6(h — h')), 
and U is the potential due to the substrate. K is the interfacial stiffness. If the wall is 
attractive, but impenetrable, U(y) = oo for y < 0, and a minimum of U at some y > 
allows to pin the interface. With long range forces, U{y) ~ ujy a ~ x +v/y a (v > 0), for large 

y !• 

The disorder due to h w and V in eq.(2) makes the partition Z[h w , V] = J Vhexp(—f3H) 



a stochastic variable. Thus, we introduce replicas |0| and evaluate 



Z n = JwJ Vh w P v [V}P w [h w ] J n Vh a exp(-{3H[h w h w ,V}) (3) 

where lia(P v ) = const — ^ / dxdhV (x , h) 2 . Integration over W is easily performed and 
allows to interpret h a (x) as world lines (x corresponding to time) of n quantum particles 
interacting via attractive two-body 5-potentials. Thus, in eq.(3) we are left with integrations 
over T>h w and T>h ai and an effective Hamiltonian: 

f3H[ha, h w ] = Jdx (^{^h w /dx p f + ^[^-(dhjdx) 2 + ^-(dh w /dx) 2 + (4) 
C(dh a /dx)(dh w /dx) + U(h a )} + AJ2 S(K - h p ) 



The logarithm of P w is now included into f3H and the couplings K' and C arise from the 
replacement h a —>■ h a + h w (initially, of course, K' — C — K). A functional renormalization 
group (RG) treatment can be performed exactly up to first order in U and A. By 
summing up exp(—/3H) over Fourier modes h a (k) and h w (k) with A/b < k < A, after the 
rescalings x — > bx, h a — > b^h a and h w — > b^ w h W) one obtains the following RG flow equations 
(b = l + dl): 

dki(U) 1 U' TJ" 

-^ = 1 + Ch u +Q li (5) 

dln(iir) 



dl 

dki(K') 

dl 
dln(C) 
dl 

d\n(K w ) 

dl 
dln(A) 



2C — 1 
2^-1 

C + - 1 

1 - 2/3 + 2C„, = 

i-C 



dl 

where Q is a suitable function of K, K w , C and the cut-off A |7]], and h stands for a generic 
h a . 

We first discuss A = 0. With ( w < 1/2 and ordered bulk, for ( = ( = 1/2 
there exists a fixed point (FP) of eqs.(5), with respect to which K' and C are irrelevant 
(d\n(K)/dl = 0, while, e.g., dhx(K')/dl < 0). Thus, substrate fluctuations decouple from 
the problem. With short range, or with long range forces such that Co — 1 /2 > 2 / \a + 1) = (* 
(strong fluctuation (SF) regime 0), the FP behavior of U at large h is within the control of 
a first order cumulant expansion and turns out to be U oc h^^, (t(Cq) = 2(1 - Co)/Co = 2) 
PJ7|JT^]. This long distance behavior should in fact apply to all the FP IPs necessary to 
describe the wetting transition in such conditions. These FP's are in general three: one 
describing pinned interface situations, one for the wet regime with unbound interface, and 
one, unstable, at the borderline between the domains of attraction of the previous two, 
describing the transition point behavior. In view of the decoupling of substrate fluctuations, 
the wetting transition controlled by these FP's, whose IPs we can not determine at finite h, 



is expected to be continuous, with exponents identical to those valid for flat wall, which are 
known exactly p|,p^|. In the case of short range forces, numerical evidence of second order 
wetting with such exponents has been recently obtained for low enough ( w by extensive 
transfer matrix calculations ||. 

The FP's for ( w > 1/2 have to be found at T = 0, by putting ( = ( w in eqs.(5). Indeed, 
now, choosing again ( = 1/2, parameters like K' and C would grow to infinity while K 
remains fixed. Surface roughness is clearly relevant now. Under a rescaling b, a T = fixed 
point is approached as (3H oc b y (j3H)* with (/3H)* finite and y > 0, when b — > oo. Such FP's 
are expected in situations when quenched disorder (due to the wall here) controls the physics 
0. At the T = FP's with ( = ( w , K, K' and C are all growing to infinity at the same rate 
( e.g., K(l) ~ K*exp[(2( w - 1)1]), and U{1) ~ U*(h) exp[(2C„, - 1)/]. U* obeys an equation 
like the first of eqs. (5), with the constant term replaced by 2(1 — ( w ), and ( w multiplying 
the second term on the r.h.s. in place of (. Thus, the discussion of the asymptotic behavior 
of U* follows lines similar to those for U in the case £ w < 1/2 0. In particular, with short 
range forces or in SF regime, we get now U*(h) oc h- T ^'\ (t(( w ) < 2). Such behavior of U* 
holds also in MF regime (C* > Co = V 2 0), as soon as ( w > (*■ 

This asymptotic behavior of U* and the connection between path integral and quantum 
mechanics are the key to demonstrate first-order wetting. Indeed, the transition order is 
revealed by the way in which (h) diverges to infinity. Consistently with eqs. (5), upon 
approaching a T = FP with ( w > ( , or ( w > (* in MF regime, we must define K w * such 
that K w (l) = K w * exp[2( w — 1)1} = const. Thus, in the FP action {(3H)* we are left with 

= 0, as / — > oo. By shifting back integration variables in this action (h a — > h a — h w ), the 
terms in K'* and C* disappear and the calculation of Z n can be easily converted into that of 
the ground state energy of a quantum problem in ID, with n + 1 particles and Hamiltonian 
H = J2 a [Pa/ '(2K*) + U(h a — h w )]. In this problem the particle with coordinate h w has an 
infinite mass. This circumstance allows to get the ground state wave function of 7i exactly 
in the form Yl a ^(h a — h w ), with \l/ satisfying the one-particle Schrodinger equation: 
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1 d 2 ^/dh 2 + U*V = ett. (6) 



(/i a — h w ) is proportional to the expectation value, (h)y W^, of h in the ground-state ^(/t) 
of eq.(0). We concluded above that, at large h and for ( w > 1/2 (or ( w > (* > 1/2 with 
long range forces in MF regime), the possible U*(h), however behaving at finite h, are 
repulsive and decay asymptotically to zero with a power t(( w ) < 2 of h. The FP U* at 
the wetting transition must have such a shape to belong to the class of borderline between 
potentials with bound ground state and e < 0, and potentials for which all states have 
e > and (h)^ = oo. These two latter types of potentials characterize dry and wet regimes, 
respectively. Independent of the details of U*{h) at short h, a solution of eq.@ with e = 
has a remarkable property for r < 2 ||10|| . Indeed, an e = eigenstate necessarily behaves as 
^(h) oc exp(— ah s ), at large h, with s — 1 — | > 0. This means that, for e = 0, a repulsive 
potential decaying slower than h~ 2 creates a too strong barrier at large distances, to allow 
interface derealization. Thus, the ground state \l/ for U* representing the transition FP 
(i.e. a FP potential in the borderline class) must be bound, with (h)^ < oo. This implies 
that, right at the wetting transition, (h) < oo, while (h) = oo as soon as the wet phase is 
accessed. First-order wetting is thus proved as soon as ( w > 1/2 ( short range forces or 
C* < 1/2), or C > C > 1/2- 

A recent numerical study of a 2D model with rough substrate exerting short range forces, 
gave evidence in support of first-order wetting for ( w sufficiently larger than 1/2 ||. In 
order to get a more direct manifestation of the mechanisms implied by eq.(6), we performed 
transfer matrix calculations for a model on square lattice with both wall and interface 
represented by directed paths, as described in ref. [nj. Fig.l reports numerical results for 



the probability distribution of h. Data are taken just below the depinning temperature 
for ( w = 2/3 (r(2/3) = 1). The dotted curve has a behavior oc exp(— ax 1 / 2 ), of the form 
expected right at threshold on an infinite asymptotic range (s = 1/2). A relatively still poor 
sampling over disorder is largely responsible of some oscillations of the distribution, but the 
overall trend appears already consistent with our theoretical predictions. 
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With bulk disorder (A > 0), the perturbative character of eqs.(5) prevents an exact 
control of the FP U for h — *■ oo. On the other hand, we know that, with A>0, Co — 2/3 is 
the exact interface anisotropy index ]]Jf3j|. For ( w < 2/3, by putting C = 2/3 in eqs.(5), we 
find that both K(l) and A(Z) grow proportional to exp(//3) (towards a T = fixed point), 
while K' and C grow slower, and are thus irrelevant. At the same time, for short range 
forces, U(l) = U* exp(//3) gives U*(h) oc h~ T ^°\ with r(Co) = 1, at large h. Thus, in the 
limit of very small bulk disorder, we get indication that for ( w < 2/3 a wetting transition 
regime identified by C = Co = 2/3 should imply a decoupling of substrate fluctuations from 
the problem. At least with very weak bulk disorder, the wetting transition with ( w < 2/3 
should retain the features of the flat wall case. For this case one indeed expects an effective 
wall-interface potential decaying as h^ 1 and Kardar has determined exactly by Bethe 
ansatz the second order character and the exponents of the transition ||. Consistently 



with our expectation, numerical results for short range forces in ref. |L6j support continuous 
wetting in Kardar's class for Q w sufficiently lower than 2/3, even with finite bulk disorder. 

Let us consider now ( w > 2/3, and short range forces again. By putting ( = ( w in eqs. 
(5), we find A(0 = A(0)exp[(l - ( w )l\, while K{1) = K'{1) = C{1) = K* exp[(2C - 1)/] 
and K w (l) = const. Furthermore, U(l) = U* exp[(2( w — 1)1} implies U*(h) oc h~ T ^ w \ Since 
now A (still supposed small) grows slower than K, C and K', it is natural to regard it as 
an irrelevant parameter with respect to the T = FP's which would be reached for A = 
strictly. Upon varying ( w > 2/3, these FP's span a subset of those already discussed with 
ordered bulk, for which quantum mechanics implies first-order wetting. Thus, we conclude 
that for C«, > 2/3 a small amount of bulk disorder is irrelevant and leaves the transition under 
control of the same mechanism outlined for pure bulk and the same ( w . Numerical results in 



ref. |I6| support this conclusion, giving evidence of first-order wetting for sufficiently large 
C w and finite disorder. Similar arguments apply to long range forces in SF and, for ( w > (*, 
in MF regime. 

In summary, our RG picture demonstrates first-order wetting in 2D with sufficiently 
rough substrates exerting short range forces on the interface. This is consistent with earlier 



numerical work suggestive of discontinuous transitions 0,0- The threshold for first-order 
wetting is precisely identified as ( w = 1/2 in the case of ordered bulk. For disordered bulk 
perturbative arguments suggest first-order as soon as ( w > 2/3, consistent with a possible 
general rule that £o identifies the threshold. We predict roughness induced first-order wetting 
also with long range forces, for ( w > (o > (* (SF) or for ( w > (* > ( Q (MF). Discontinuous 
depinning is due to the repulsive effective wall-interface potential, which becomes too strong, 
at large distance, to allow for a continuous increase towards infinity of (h) when depinning 
is approached. This follows from general quantum properties, independent of the details of 
U at finite h. 

Interesting open problems remain the nature of wetting right at the thresholds and the 
possible extension to 3D of this type of results, which rely on the connection with quantum 
mechanics in ID. A recent mean feld study in 3D suggests the possibility of first-order 
wetting induced by wall roughness with short range substrate potential and ordered bulk 
|TT|| . Another interesting issue is whether ( w = ( could be a plausible threshold also 
in cases in which different kinds of bulk disorder imply different £o' s - Relevant examples 
include random fields [|j] and quasicristals [|T^ . 

Due to the competition between two qualitatively similar scaling geometries, interactions 
between a fluctuating manifold and a random boundary can lead to interesting phenomena 
also in other contexts. An example could be flux lines in high-T c superconductors with 
extended rough defects [19j . Also of interest would be polymers or membranes adsorbed by 
rough walls. 
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FIGURES 

FIG. 1. Inset: sketch of the geometry of wall (continuous) and interface (dashed). Main: 
probability distribution of h, from a sampling of 500 wall configurations of 100000 longitudinal 
steps. 
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